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МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ 
РАЗВЕТВЛЕННЫХ ГИДРАВЛИЧЕСКИХ СЕТЕЙ 


Рассмотрена математическая модель сложных разветвленных гидравлических Сетей, позволяющая рассчи- 
тывать распределение жидкости во всех точках. 
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Введение. Водораспределительные сети городов — сложные системы, включающие насосные 
станции, водонапорные башни, насосы подкачки и трубопроводы, состоящие из множества труб 
различного диаметра и длины и большого числа потребителей. Возможности современной вы- 
числительной техники позволяют проводить анализ работы сетей любой сложности. 
Постановка задачи. Некоторые авторы расчет распределения жидкости в разветвленной се- 
ти предлагают проводить по упрощенной динамической модели без учета инерционных сил 
жидкости и волновых процессов. Такая модель не дает правильной картины переходных про- 
цессов, но при отсутствии кавитационных явлений на участках сети обеспечивает достоверный 
результат в установившемся режиме работы. 

Основой программы модели сети являются математические модели всех входящих в нее 
элементов. Рассмотрим работу одного участка сети (рис.1). 
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Рис.1. Расчетная схема участка сети (В ,4\!›2! - давление, расход и геометрическая высота одного конца линии 1 
(первый индекс), входящего в узел 1 (второй индекс); В, 412,2.) - давление, расход и геометрическая высота второ- 


го конца линии 1, входящего в узел 2; [| ‚4,5, - длина, внешний диаметр, толщина стенки трубы участка 1) 


Решение задачи. Условимся считать расход, входящий в линию, положительным, а вытекаю- 
щий — отрицательным. 

Без учета инерционных сил жидкости и волновых процессов работу одного участка мож- 
но описать следующими уравнениями: 





Ч =/(Р,); (1) 

4 = шо | Авар ; (2) 
В, +р2, =Р,› +р82, +АР; (3) 
ИС (чи +4). (4) 
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Уравнение (1) - это характеристика источника (насоса). Уравнение (2) представляет со- 
бой характеристику стока (потребителя), где и — коэффициент расхода, « — площадь отвер- 
стия, р — плотность жидкости, д — ускорение свободного падения. 

Уравнение Бернулли (3) записано для узлов 1 и 2, где аР; - потери давления в линии 1, 
при этом 

Г ру? 
АР = р 
а 2 


где / — длина трубы; /, - коэффициент гидравлического трения; 4 - внутренний диаметр тру- 
бопровода; у — средняя скорость жидкости в трубе. 
Уравнение (4) связывает производную давления по времени с приведенной жесткостью 
С, рассматриваемого участка и балансом расходов по концам линии, при этом 
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где Ет› — модуль упругости материала трубопровода; Е» -— модуль упругости жидкости; 
4, - наружный диаметр трубопровода; 5’- толщина стенки трубы. 
Таким образом, система уравнений (1)-(4) позволяет определить параметры потока на 
каждом участке. 


Основываясь на этих уравнениях, можно составить математическую модель сети любой 
сложности. Рассмотрим сеть, состоящую из двух насосов, двух потребителей и четырех линий 


(рис.2). 
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Рис.2. Схема двухнасосной сети 


Расходы а, на стороне линии 1, присоединенной к узлу 1, и 45 на стороне линии 2, 


подсоединенной к узлу 2, определяются характеристикой источников и зависят от давлений в 
узлах: 


4и = (В); (5) 
4»› = Л(Р,). (6) 


Расходы а.., на стороне линии 3, подсоединенной к узлу 3, и 4. на стороне линии 4, 
подсоединенной к узлу 4, соответственно: 


|2 | 

Чзз = Из РТР) р (7) 
2 . 

444 = Иа Риме Р) . (8) 
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Уравнения Бернулли для линий: 








В, +р82, = ВА; +р825 +АВ; (9) 
Р,› +р82, = Р; +р82; + 4Р,; (10) 
Р; +ря2; = В; +р82. + АР; (11) 
Р,; +ря2; = Ри +ре2. + АР. (12) 
Производные давления в линиях: 
АР 
а. +415); (13) 
АР 
— = СЧ» +45); (14) 
АР. 
—- = С; (435 + 433); (15) 
АР 
— = С. (445 +44). (16) 
Баланс расходов в узле 5: 
415 + 455 + 435 +445 =0. (17) 


Следствием соединения линий в узле 5 является равенство давлений на сторонах линий, 
подсоединенных к узлу, давлению в узле 5: 

В. =Р;=В; =Р,;=Р. (18, 19, 20, 21) 

Таким образам, мы получили систему из 17 уравнений, которые определяют 17 неиз- 

вестных. Четыре уравнения дифференциальные (число линий), остальные алгебраические. Рас- 

ходы источников (411,4) и стоков (4..,4..) выражены в явном виде. Получим выражение в 





явном виде для расходов на сторонах линий, сходящихся в узле 5 (415,45, 935, Ч4 ). 


Продифференцируем выражения (9)-(11) по времени и, пренебрегая величинами второй 
степени малости (при решении системы уравнений численным методом изменение расхода на 
шаге интегрирования есть величина бесконечно малая), имеем: 








АР =: АР; ь (22) 
@& Ш’ 

аР.› Ее аР,; : (23) 
@& @' 

аР;; ра АР; р (24) 
[ИА [ИА 


Заменяя в уравнениях (18)-(20) давления на их производные с учетом уравнений 
(22)-(24) получим три уравнения, связывающие расходы на сторонах линий, сходящихся 
в узел 5: 


С! (Чи +415) = С, (4»› +94); (25) 
С› (4.2 +4») = С. (435 +4зз); (26) 
С, (435 +4зз) = С. (445 +4). (27) 
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Решая систему алгебраических уравнений (25)-(27), (17) находим в явном виде расходы 
на сторонах линий в узле 5: 





1 

Ч15 = (Ч +922 + 433 +4и)-Чи; (28) 
1 1 1 ГС 

п | 


ее. 





1 
45 = ] ] ] ГС, (Чи +422 + 433 +4и)-4>; (29) 











4з5 = ] ] ] Г@ (аи + 422 1 433 +4и)—4зз й (30) 





(Ч +422 + 433 +4и)-4щ. (31) 





945 = 
1 й 1 1 й С, 
С С, С) С. 
Следовательно, работа сети (рис.2) описывается системой уравнений (5)-(16), (21), 
(28)-(31), решение которых численным методом позволяет определить все неизвестные. 
Можно записать общую формулу определения расходов на стороне и линий, подклю- 
ченных к А-му узлу: 





1 1 

Ч ль = не м. (32) 
пс, 

С 


=] 


Здесь 4 „‚) - расход на стороне /-ой линии, подключенной к А-му узлу; С’, - жесткость /-ой ли- 
нии; 4х.) - Расход на другой стороне /-ой линии. 


Заключение. Таким образом, уравнения (1)-(4), (32) являются основой для математической 
модели водораспределительной сети любой сложности и различной топологии. На основе мате- 
матической модели составляется программа решения системы уравнений численным методом с 
использованием подпрограммы Рунге - Кута, которая позволяет анализировать на ЭВМ работу 
водораспределительных сетей любой сложности, но в связи с принятыми допущениями не по- 
зволяет анализировать переходные процессы в сети. Для их анализа необходима уточненная 
математическая модель, учитывающая инерционные силы жидкости, волновые и кавитационные 
явления. 


Материал поступил в редакцию 11.11.10. 
\.$. КВОТТКОХ, К.А. Г1КНОБЕРО, М.\. ОТУЦЫМ 


МАТНЕМАТТСАЕ МОБЕ ОЕ 5ОВОТУТОЕО НУОВАЧЕТС СТВСИТТ$ 


Тре таетайса/ тодЕ!/ оЁ зибамаеа вуагаи!с сгсийз паЁ регтйз сасшайоп о! пе Пиа @тБивоп т а! 
Роб 15 сОП$!АЕГеа, 
Кеу игога&: Вудгаи!с сгсий5, те{тетайса!/ тодЕ/та, сасшайоп о! пи @тБийоп. 
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